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Recent developments of microscopic mechanical experiments allow the manipulation of individual polymer 
molecules in two main ways: uniform stretching by external forces and non-uniform stretching by external 
fields. Many results can be thereby obtained for specific kinds of polymers and specific geometries. In this 
work we describe the non-uniform stretching of a single, non-branched polymer molecule by an external field 
(e.g. fluid in uniform motion, or uniform electric field) by a universal physical framework which leads to 
general conclusions on different types of polymers. We derive analytical results both for the freely-jointed 
chain and the worm-like chain models based on classical statistical mechanics. Moreover, we provide a Monte 
Carlo numerical analysis of the mechanical properties of flexible and semi-flexible polymers anchored at one 
end. The simulations confirm the analytical achievements, and moreover allow to study the situations where 
the theory can not provide explicit and useful results. In all cases we evaluate the average conformation of 
the polymer and its fluctuation statistics as a function of the chain length, bending rigidity and field strength. 



I. INTRODUCTION 

Modern methods for stretching single molecules pro- 
vide a valuable insight about the response of polymers 
to external forces. The interest on single molecules load- 
ing encouraged new research and technological develop- 
ments on related mechanical experiments. Typically, me- 
chanical methods allow the manipulation of a polymer 
molecule in two ways: the stretching of the chain by the 
direct action of an external force or by the application of 
an external field. If we consider homogeneous polymers 
(with all monomers described by the same effective elas- 
tic stiffness), then we obtain a uniform strain with the 
external force and a non-uniform strain with the applied 
field. 

To exert an external force on a polymer fixed at one 
end, laser optical tweezers (LOTs)^, magnetic tweezers 
(MTs)^ or atomic force microscope (AFM)^ can be used. 
Many experiments have been performed over a wide class 
of polymers with biological relevance, such as the nucleic 
acids (DNA, RNA)**, allowing the stretching of the en- 
tire molecule and providing the reading and the mapping 
of genetic information along the chaini^ Furthermore, 
it has been possible to describe the elastic behaviour of 
single polymers consisting of domains which may exhibit 
transitions between different stable stateSf^"— Other in- 
vestigations performed on double-stranded DNA deter- 
mined the extension of the polymer as a function of the 
applied forceiS, providing results in very good agreement 
with the Worm-Like Chain (WLC) modeUi"— and the 
Freely- Jointed Chain (FJC) modelJ^^ii 
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Alternatively, it is possible to manipulate single 
molecules by an external field. In this case the external 
field acts on the molecules from a distance or, in other 
words, without a defined contact point for applying the 
traction. A non-uniform stretching performed by an ex- 
ternal field can be induced either via a hydrodynamic (or 
electrohydrodynamic) flow fieldi^^— or via an electric (or 
magnetic) fieldJ^"— One experimental advantage of us- 
ing flow fields is that the liquid surrounding the tethered 
molecule can be easily replaced; this is indeed an im- 
portant feature for many single-molecules studies of en- 
zymes which require varying buffer conditions.— The flow 
field technique was extensively applied in single-molecule 
study of DNA elasticity^ ^ as well as to characterize the 
rheological properties of individual DNA molecules<^"— 
The use of an electric field has been adopted for driving 
the alignment of DNA on a solid surface for applications 
such as gene mapping and restriction analysis. ■'^^ Finally, 
magnetic fields have been used to apply torsional stress 
to individual DNA molecules 

In order to understand the response of polymers to ex- 
ternal fields and to study their statistics, some theoretical 
models have been proposed. These models are typically 
based on the FJC and WLC schemes, generalized with 
the inclusion of the given applied field. Some studies have 
shown that in a weak external field the persistence length 
along the field direction is increased, while it is decreased 
in the perpendicular direction; moreover, as the external 
field becomes stronger, the effective persistence length 
grows exponentially with the field strength.—"— Other 
investigations under a constant velocity flow have shown 
that a flexible polymer displays three types of confor- 
mation: unperturbed at low velocity; "trumpet" shaped 
when partially stretched; "stem and flowers" shaped. 
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with a completely stretched portion (the stem) and a se- 
ries of blobs (the flowers) , at larger loading^"^ Polymer 
models have been studied in elongational flows to analyze 
the coil stretching and chain retraction as a function of 
polymer and flow parameters, finding good agreement 
with experimental datai^i^ Conformational properties 
of semiflexible polymer chains in uniform force field were 
also studied for two-dimensional models."^'' In spite of all 
these relevant efforts, it is yet a challenge to base on one 
same unified theoretical framework and understanding of 
all aspects of polymer mechanics in an external field. 



Building on our previous studies jSiii^ in this paper we 
study the conformational and mechanical properties of 
flexible and semi-flexible non-branched polymer model 
chains tethered at one end and immersed in an external 
force field. This situation is useful to describe almost two 
physical conditions of interest: a polymer chain immersed 
in a fiuid in a uniform motion (our model is valid only 
when the action of the fiuid motion can be described by 
a distribution of given forces applied to all monomers) 
and an arbitrarily charged chain inserted in a uniform 
electric field. 



Our theoretical approach is twofold, since we adopt 
both analytical (statistical mechanics^^i^) and numer- 
ical techniques (Monte Carlo simulations^SiS) . While 
the analytical approach is useful to obtain the explicit 
partition function in some specific cases, Monte Carlo 
simulations are crucial to study more generic cases, in- 
accessible to analytical treatments. In particular, while 
we develop our theoretical framework starting from the 
more tractable FJC model, we take full profit from our 
MC simulations to extend our study also to the WLC 
model. 



The structure of the paper is the following. In Section 
II we introduce the mathematical formalism adopted and 
we derive a generic form of the partition function in 
for a generalized FJC model where the extensibility of 
the bonds is taken into account. In the Section III we 
find the two specific forms of the partition function for 
the 2D- and the 3D-case for the pure FJC polymer with 
non extensible bonds. Moreover, we obtain in both cases 
the variance and the covariance among the positions of 
the monomers. In the Section IV we present the gen- 
eralization of previous results to the semi-flexible WLC 
model. We present two closed-forms approximations for 
the 2D- and the 3D-case and the comparisons with MC 
simulations. In section V we analyze the behavior of a 
chain in an external fleld to which also an external force 
is applied at the end of the chain. The case with the force 
not aligned with the fleld is particularly interesting and 
shows the power of the MC method. Finally, in Section 
VI some conclusions are drawn. 
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FIG. 1. (color online) A polymer chain in an external field. 
The first monomer is clamped at position rb while the oth- 
ers are free to fluctuate. Each monomer is subjected to an 
external force qk (different in strength and direction for any 
K): all these forces mimic an external field. Another external 
force, playing the role of a main pulling load, /, is applied to 
the last monomer at the position rjv. 



II. GENERAL THEORETICAL 
FRAMEWORK 

As previously discussed, the polymer models most used 
in literature are the FJC and the WLC. As argued in 
Rcf,38, for weak tension and weak external field, it is ac- 
ceptable to model the polymer as a FJC model. This 
model breaks down only when the curvature of the con- 
formation is very large because it ignores the consequent 
great bending energy. Since we will look upon this prob- 
lem in the end of this work, we now give way to the 
case of a FJC. In particular we consider a FJC with two 
additional hypothesis. Firstly we consider the possible 
extensibility of the bonds of the chain through a stan- 
dard quadratic potential characterized by a given equi- 
librium length: such an extension mimics the possible 
stretching of the chemical bond between two adjacent 
monomers. If necessary, the extensibility of the bonds, 
here described by linear springs, can be easily extended 
to more complex, nonlinear springs. "^^ Moreover, we take 
into account a series of arbitrary forces applied to each 
monomer: these actions mimic the effects of an external 
physical field applied to the system. In addition, we con- 
template the presence of an arbitrary force applied to the 
terminal monomer of the chain. All calculations will be 
performed in 3?'' and we will specialize the results both 
in the 2D-case and in the 3D-case when needed. The 
idea is to write the complete form of the Hamiltonian of 
the system and to build up the corresponding statistical 
mechanics. ^■^ The starting point is therefore the calcula- 
tion of the classical partition function. In fact, when this 
quantity is determined, it is possible to obtain the force- 
extension curve (the equation of state) through simple 
derivations. 

Let us consider a non-branched linear polymer with N 
monomers (see Fig. [T|) at positions defined by ri, r/v G 
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5ft'' (for considering d = 2 or d = 3 according to the 
specific problem of interest). To eaclr monomer a given 
external force is applied and named gi,...,gN. Another 
external force, playing the role of main pulling load, /, is 
applied to the last monomer at the position f^- While 
the chain is clamped at position fo, the monomers are free 
to fluctuate. The Hamiltonian of the system is therefore 
given by 



N ^ ^ N 



N 



where pi are the linear momenta, m the mass of the 
monomers, k the spring constant of the inter-monomer 
interaction, and I the equilibrium length of the monomer- 
monomer bond. We search for the partition function of 
the system defined as: 



exp 



H 



dri . . . dr jv dpi ■ • ■ dpW ( 2 ) 



2Af-times 



where c is a multiplicative constant which takes into ac- 
count the number of microstates. As well known, the 
kinetic part can be straightforwardly integrated and it 
yields a further non-infiuencing multiplicative constant; 
then we can write the partition function as an integral 
over the positional space only. This integral can be easily 
handled through the standard change of variable 



' ^1 = f 1 - fo 

£,2 = r2- n 

. = rN - rN- 



(3) 



having the Jacobian determinant J 



1. 



3(^1. ..fjv) 
c?(Ci...«Jv) _ 

We consider the terminal fg of the chain fixed in the 
origin of axes, i.e. fo = 0. So, we cast the positions f in 
terms of the variables ^ as follows 



ri 



n=^2 



(4) 



By setting the general solution as f; = X]k=i^^' 



partition function becomes 



Zd^ c I ...I exp 



N 



jv-times 
X exp 



2kBT 



K 
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K=l 



K 



(5) 
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X exp 
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1 

d^i...d^N 



K=l ,J=1 
N 



^/ • E 

Inverting the two summation symbols 

N K N N 

E • E ^"-^ = E ■ E 9^ 



we obtain 



where 



K=l 



Zd — 



J=i 



K=l 



i=K 



N „ 



2kBT 



> 



(6) 

(7) 

(8) 
(9) 



It exists a deep conceptual connection between the last 
integral for the partition function and the theory of the 
d-dimensional Fourier transforms. The Fourier integral 
of an arbitrary function /(f) is defined as 



F(^) = / /(Oe-^-«de 



with inverse transform given by 
If we consider 



(10) 



(11) 



(12) 



it is easy to realize that the integral in Eq.(I7|) is the 
Fourier transform of /(f) calculated for uj — iVk: i-G- 



N 



Zd^c\[ F{iVk) 



(13) 



K=l 



with a e Vk defined respectively in Eq.® and Eq.®. It 
is important to remark that the function in Ea. (jl2p has 
a spherical symmetry (i.e. it depends only on the length 
of the vector £) and, therefore, also its Fourier transform 
F[u) exhibits the spherical symmetry, depending only 
on the quantity |aJ| in the transformed domain. In fact. 
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for such spherically-symmetric functions it holds that: if 
/(^) = /(l^) then = F{\uj\). Furthermore, we have 
that 

F{n) = £^27rpf{p)(^y J._,{pn)dp (14) 



for d = 2n (even), and 



+00 



\ n J 



/2TTPY 



Fi^)=J^ 4V/(p)(^^j Ji_l{p^)dp (15) 

for d = 2n + 1 (odd), where p = \^\ and = 
Here J,^{z) and jiy{z) are the cylindrical and spherical 
Bessel functions of the first kind respectively, correlated 
by the standard relation jtx(z) — y^^J^^i (z)4ii^In our 

calculations we have to set lo = iVk and, therefore, we 
obtain Q, — i\VK\- Moreover, when the argument of J,y(2:) 
and ji/(z) is supposed imaginary we obtain the modified 
Bessel functions of the first kind^i^ 



i^{z) = {iyj^{iz) 



(16) 



For example we have the explicit expression jo(^) — 
and io{z) = sii^ while, on the contrary, Io{z) and Jq{z) 
cannot be written in closed form. So, for d even we even- 
tually obtain 



III. FREELY-JOINTED CHAIN MODEL 
UNDER EXTERNAL FIELD 

A. Average values of positions 

In the previous section we obtained the general ex- 
pression of the partition function for the case where the 
extensibility of the bonds is taken into account. This is 
described by the parameter k, which characterizes the 
elastic bond between adjacent monomers. In the present 
Section we want to study the effects of an arbitrary dis- 
tribution of forces on a pure freely jointed chain model 
(FJC). Therefore we need to obtain the specific form 
of the partition function in the case of rigid bonds of 
fixed length I. From the mathematical point of view it 
means that we will consider /c — > cx), a condition rep- 
resenting a inextensible spring. Because of the relation 
-ax _ ^^^-j -y^fjjgQ a — >■ 00 we may determine the 



limit of Eg. dTOl) and Eg. ((20)) for a — >• 00 (i.e. for fe — )■ 00, 
FJC limit). Since the arbitrariness of the constant c, 
we may consider in Eqs. and ([^0]) a multiplicative 

constant term (\/^)^- Then, by using the translated 

property -y/f e~°^''~'-' 5{p — /) for a — >■ 00 we perform 
all the integrals thereby obtaining 

^ 1 

Zd = cW ^ ^ Ii^{l\VK\) deveii (21) 



FjiVK) 
27r 



\Vk\ 



and, on the other hand, for d odd we have 



(17) 



47r 



id^{p\VK\)dp 

Finally, by using Eg. (1131) . the partition function is given 
by 



Id^{p\VK\)dp 
(19) 




for d even, and 



N 
K=l 



+00 



(20) 

for d odd, where a and Vk are given in Eqs.® and ([9]). 
In the framework of statistical mechanics, the knowledge 
of the partition function allows to determine all needed 
expected values describing the statistics of the chain (i.e., 
average values of the positions, variances of the positions 
and so on). 
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In particular, for d = 2 we have 



d odd 



N 



K=l 

while for d = 3 wc obtain 



K=l 



N 



/ + E. 



sinh 


[kBT 


f + J2f=K 9i 


) 


k 


I 

bT 


f + J2^=K 9^ 





(22) 



(23) 



(24) 



All the expressions given in Eqs.dH]), (gl]), (ESI, (HH) can 
be summarized in the general form 



N 



Zd = cl[fi\VK\) 



(25) 
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with a suitable function f{x). By using this expression of 
the partition function wc can find the average position of 
the z-th monomer of the chain; indeed, from the definition 
of the Hamiltonian in Eq.(IT]) we state that fi — — and, 
therefore, we get 



d 

in) = ksT— \nZd 
09% 



(26) 
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which represents the shape of the polymer chain under 
the effects of the external field cji and the applied force /. 
Now we can substitute Ea.ipSj) into Eg. ipS]) . obtaining 



E 



Vk 
\ \Vk\ 



1 df{x) 



f{x) dx 



(27) 



--\Vk\ 



In 2D we have f{x) = Iq{Ix) and therefore we obtain 



K=l Iq 



I 



f 



2^j=k9j 



f + Y.J=K 9 J 



f + J2j=K 9 J 

(28) 

For such a 2D case, by applying Ea. ((28|) . the average val- 
ues of the longitudinal component of the positions have 
been calculated and are plotted in FigH] as a function of 
the chain length N and the field strength g. We have 
considered only the action of an external uniform field 
with gj — g and amplitude g. 

Although this case lends itself to a full analytical so- 
lution, numerical simulations were also performed by us- 
ing a conventional implementation of the Metropolis ver- 
sion of the Monte Carlo algorithm.''^ The initial state of 
the chain is defined by a set of randomly chosen posi- 
tions. The displacement extent of each step governs the 
efficiency of the configurational space sampling. There- 
fore, we analysed several runs in order to optimize its 
valuei^i^ The perfect agreement between the theory and 
the MC simulations provides a strict check of the numer- 
ical procedure, to be used in the foregoing. 

On the other hand, in 3D we have f[x) =- 
leading to 



Ix 



I 



K=l 



N 



g.j 



f + Ej=K 9.1 



2^,7= 



k9.i 



(29) 

where £{x) = cotha; — ^ is the Langevin function. By 
using Eq. ip^ . as before, it is possible to plot the aver- 
age values of the longitudinal component of the positions 
for the 3D case (Fig[2|). Also in this case we adopted 
a uniform field g and the good agreement with the MC 
simulations is evident. 

As particular case, if there is only the force / applied to 
the system we obtain the standard scalar force-extension 
curves linking r — |(rw)| with / = |/|. In 2D we have 



r 

In 



\kBT J 



(30) 



in agreement with recent results^^^ while in 3D we obtain 



IN 



If 



(31) 



which is a classical result.— The simple results in 
Eas. ([5ni) and (I5T]) have been used to obtain the limiting 




FIG. 2. (color online) Average values of the longitudinal com- 
ponent of the positions induced by the external field for the 
2D FJC case. The red solid lines correspond to the analyt- 
ical results Eas. H28|l and ([32}, MC results are superimposed 
in black circles. Top panel: each curve corresponds to dif- 
ferent chain lengths N = 10, 20, 30, 40, 50 for a fixed value 
gl/(kBT) = 1 (e.g., corresponding to / = Inm, g — 4pN at 
T = 293K). Bottom panel: each curve corresponds to the 
different values gl/iksT) = 0.1,0.25,0.5,1,2,10 for a fixed 
chain length iV = 20. 



behaviors under low (/ 0) and high (/ — > oo) values 
of the applied force, as shown in Table ID 

Building on such first results we now focus on some 
particular interesting approximations. More specifically, 
it can be interesting to find approximate results for the 
case of a homogeneous field and no end-force, f — and 
gj = g for any J. In this case we search for the scalar 
relation between r = [("nv)! and g = \g\. In the 2D case, 
from Eg. (1^51) . we have 



^^^^/i(^(A^-fc + l)) 
IN N 



1 

N 



fe=i/o(^(A^-fc + l)) 
^Ii{^{N-x + 
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NJsI 



° ^o(^(iV-x + l)) 



dx 



log- 



(32) 
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FIG. 3. (color online) Average values of the longitudinal com- 
ponent of the positions induced by the external field for the 
3D FJC case. The red solid lines correspond to the analyt- 
ical results Eas. H29|l and ()33|l . MC results are superimposed 
in black circles. Top panel: each curve corresponds to dif- 
ferent chain lengths N = 10, 20, 30, 40, 50 for a fixed value 
gl/iksT) = 1. Bottom panel: each curve corresponds to the 
different values gl/iksT) = 0.1,0.25,0.5,1,2,10 for a fixed 
chain length iV = 20. 



B. Covariances and variances of positions 

In this Section, we search for the covariance among 
the positions of the monomers. It is important to evalu- 
ate such a quantity in order to estimate the variance of 
a given position (measuring the width of the probabil- 
ity density around its average value) and the correlation 
among different monomer positions (measuring the per- 
sistence of some geometrical features along the chain). 
In order to do this, we identify the a-th component of 
the i-th monomer as . The covariance of the generic 
monomer simply defined as (it represent the expectation 
value of the second order): 

Cov(ria,rj^) = {{r^a - {ria)){rji3 - (rjis))) (34) 
= inarjis) - {ria){rji3) 

Taking the derivative of the partition function with re- 
spect to the a and the /3 components of the force vectors 
Qi and gj we can solve the problem as follows. We con- 
sider the standard expression for the partition function 
and we can elaborate the following expression 



dgia dgjfi dg^adgjp 
or, equivalently, by introducing Ea.([26| 

d 



but we can simply determine that 
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x = \Vk\ 



On the other hand, for the 3D case we obtain 
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{N + 1) (e 



-1 (33) 



We have usefully exploited the fact that, for large N, the 
sums can be approximately substituted with the corre- 
sponding integrals, which are easier to be handled. The 
closed- form expressions given in Eas. (l32p and p3p are 
very useful to obtain the limiting behaviors of the poly- 
mer under low (g — > 0) and high (g — >■ oo) values of the 
applied field, as shown in Table HI Moreover, we have 
verified the validity of Eos. (15^ and ([55)1 through a series 
of comparisons with MC results (see Figl6] in the next 
Section for details). 



where we have defined the unit vector cTq as the basis of 
the orthonormal reference frame. Being 



if ^ \ 
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and 
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d\VK\ 

dgjp ksT \Vk\ 
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q=K 



Jq 
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after long but straightforward calculations we obtain 
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Ordering the terms we finally obtain the important result 



Cov{ria,rjp) 



6., f{\VK\) 
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\Vk\ fi\VK\) 



(41) 
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It represents the final form of the covariance between 
two different components of the positions of two different 
monomers. 

If we look at the variance of a single component of a 
single position (i = J, a = 13) we have the simpler result 
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f{\VK\) 
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K=l \VK\f(.\VK\) ^1 \VK\'f{\VK\) 

( f'i\VK\) fi\VK\r\ 



(42) 



X /"(lV^i.1) 



\Vk\ fi\VK\) 



In order to use the previous expressions we have to 
specify the function / and its derivatives for the two- 
dimensional and the three-dimensional case. In the 
2D case we have f{x) = Iq{Ix), f'{x) — lli{lx) and 
f"[x) = ^^[Iq{Ix) + l2{lx)]. On the other hand, for the 
3D case we have f{x) ~ ^"^^^^^ , f'{x)/f{x) ~ lC{lx) 
and f"{x)/f{x) = P - 2lC{lx)/x. This completes the 
determination of the covariance. 

We report in FigU] and FiglS] the longitudinal and 
transversal component of the variance as a function of the 
chain length and the field strength for the 3D case (with 
/ = 0). The 2D case is very similar and it has not been 
reported here for sake of brevity. We can observe some 
interesting trends: the longitudinal variance of the posi- 
tion is a decreasing function of the number of polymers 
N while the transversal one is a increasing function (with 
a fixed amplitude of the external field g). Moreover, both 
variances are rapidly increasing along the chain, assum- 
ing the largest value in the last free monomer, which is 
more subject to strong fluctuations. It interesting to ob- 
serve that the variance (both longitudinal and transversal 
components) is a linear function of the position i along 
the chain (it linearly intensifies along the chain itself) 
with a simple force / applied at the free end: conversely, 
with a uniform field the distribution of forces generates 
a strongly non-linear intensification of the variances mov- 
ing towards the free end-terminal. So, from the point of 
view of the variances, the application of a field or the ap- 
plication of a single force generates completely different 
responses. In Fig|5] we can also observe that the vari- 
ances are decreasing functions of the strength of the field 
(both for the longitudinal and transversal components); 
in fact, the intensity of the fields tends to reduce the fiuc- 
tuations of the chain, increasing, at the same time, the 
tension within the bonds. 




FIG. 4. (color online) Longitudinal (top panel) and transver- 
sal (bottom panel) component of the variance of positions for 
the 3D FJC case. The red solid lines correspond to the ana- 
lytical result Eq. (|42p . MC results are superimposed in black 
circles. Each curve corresponds to different chain lengths 
iV = 10, 20, 30, 40, 50 for a fixed value of the external field 
defined by gl/iksT) = 1. 



IV. WORM-LIKE CHAIN MODEL UNDER 
EXTERNAL FIELD 

In previous Sections we treated systems described by 
the FJC model, characterized by the complete flexibil- 
ity of the chain and, therefore, by the absence of any 
bending contribution to the total energy. Neverthe- 
less, in many polymer chains, especially of biological ori- 
gin, the specific flexibility (described by the so-called 
persistence length^) has a relevant role in several bio- 
mechanical processes. In order to take into consideration 
these important features, with relevant applications to 
bio-molecules and bio-structures, in this Section we in- 
troduce the semi-flexible polymer chain characterized by 
a given bending energy added to the previous Hamilto- 
nian 



N 



^ = E 



2m 



N 



\rK-rK-i\\-l) (43) 



K=l 



^ N-1 N 

-K {U+i ~ U) -^gK-fK-f-rN 



i=l 



where k is the bending stiffness, k is the stretching mod- 
ulus and U = (r^+i — ri)/||r'i+i — ri\\ is the unit vector 
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T{fl/{kBT)). When / = and gj = g for any J we 
search for the 2D scalar relation between r and g = \g\. 
As discussed in a previous section (see Eas.(|32)) and (|33)) ). 
we can write 



1 ^ 



k=l 



{N-k + 1) 



k=0 



ig 

knT 



iN-x + l)]dx 



1 1 



^ (y) dy 



(45) 



where we have defined the change of variable y = 
-^^{N — x + 1). We adopt now a second change of vari- 
able through the relation z = F{y) ot y = J-^^{z); it 
leads to 



r 

In 



1 1 



1 1 / 



dz 



dz 



(46) 



FIG. 5. (color online) Longitudinal (top panel) and transver- 
sal (bottom panel) component of the variance of positions for 
the 3D FJC case. The red solid lines correspond to the analyt- 
ical result Ea. (|42|l . MC results are superimposed in black cir- 
cles. Each curve corresponds to different values of the external 
field amplitude defined by gl/iksT) = O.f, 0.25, 0.5, f, 2, 10 
for a fixed chain length A'^ = 20. 



collinear with the i-th bond (see lief Hi for details). In 
particular we take into consideration the classical WLC 
model, describing an inextensible semi-flexible chain: it 
means that the spring constant k is set to a very large 
value (ideally k — ^ oo) so that the bond lengths remain 
fixed at the value I. It is well known that it is not possible 
to calculate the partition functions in closed form for the 
WLC polymers. Nevertheless, some standard approxi- 
mations exist for such cases leading to simple expres- 
sions for the force-extension curves when a single force / 
is applied to one end of the chain. In the following, start- 
ing from these results, we search for the force-extension 
curves when the polymers is stretched through a constant 
field g. 

We start with the result for the 2D- WLC with an ap- 
plied force /: the approximated force extension curve is 
given by^ 



fl 
kuT 



I 



1 _ l_ 7 ■ 
16(1 -C)2 16^8^ 



(44) 



where C, — r/{lN) is the dimensionless elongation and 
Lp = ln/{kBT) is the persistence length . We suppose 
that such a constitutive equation is invertible through 
the function leading to the expression Q = r/(lN) = 



1 



1 



16 



8(1 - z) 16(1 - z)2 



where we used the notation [^.(x)]^ = h{b) — h{a). This 
result represents (although in implicit form) the approx- 
imated force-extension curve for the 2D- WLC under ex- 
ternal fields. To evaluate Ea. (l46l) we need to know the 
inverse function J^(-), a task that can be performed nu- 
merically. 

Similarly, we may consider the standard 3D- WLC 
model with an applied force /; the classical Marko-Siggia 
resulti^ is 



fl 
kuT 



I 



1 



4(1-0^ 



1 



(47) 



where, as before, C = r/{lN) is the dimensionless elon- 
gation and Lp = lK/{kBT) is the persistence length. We 
suppose again that such constitutive equation is invert- 
ible through the function leading to the expression 
C - r/{lN) - g{fl/{kBT)). When / = and gj = g 
for any J wc search for the 3D scalar relation between 
r and g = \g\. By repeating the previous procedure, we 
can write 



1 



IN N J9- Jal '9 

1 1 



dz 



dz 



(48) 
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FIG. 6. (color online) Force-extension curves of a FJC poly- 
mer in an external field (or external force) with N=20. The 
red line corresponds to the approximated expressions given 
in Eas. (|32|l and Eas. (|33|) while the black circles have been 
obtained through MC simulations. The 2D (Eq.(0) and 
3D fEa. pi|) ) FJC expressions (without an external field) are 
plotted for comparison with f — g and / — Ng. 



which represents the impHcit form of the approximated 
force-extension curve for the 3D- WLC under external 
fields. 

It is interesting to compare the very different force- 
extension curves for a single molecule in the two cases 
of a uniform (only / applied) and non-uniform (only 
g applied) stretch. In particular, taking advantage of 
our approximated formulas, we can analyse the case of a 
FJC and a WLC polymer. The 2D and 3D FJC resuhs 
are plotted in FiglH] on the other hand, the 2D and 3D 
WLC curves have been shown in Fig[71 For the WLC 
case we assumed k = lOfc^T for the bending modulus at 
T = 293K. This value is comparable to that of polymer 
chains of biological interest (e.g., for DNA k = ISfesT).^^ 
In any case three curves have been reported for drawing 
all the possible comparisons: the response under the field 
g, the response under the force f = g and, finally, the re- 
sponse to an external force / = Ng. Interesting enough 
we note that the curve corresponding to the field g is 



FIG. 7. (color online) Force-extension curves of a WLC poly- 
mer in an external field (or external force) with N=20. The 
red line corresponds to the approximated expressions given in 
Eas. (|46|l and Eas. H48|) while the black circles have been ob- 
tained through MC simulations. The 2D (Eq.^) and 3D 
(Eq.(|47|) WLC expressions (without an external field) are 
plotted for comparison with f = g and / = Ng. The value of 
the bending spring constant is «: = 0.4 • 10~^^ Nm ~ lOksT 
at T = 293K. 



always comprised between the cases with only the force 
f — g and / = Ng. The response with the field g is 
clearly larger than that with the single force f = g since 
the field corresponds to a distribution of N forces (of in- 
tensity /) applied to all monomers; therefore, the total 
force applied is larger, generating a more intense effect. 
However, the case with a single force / = Ng shows a 
response larger than that of the field g. In this case the 
total force applied in the two cases is the same but the 
single force Nf is applied entirely to the last terminal 
monomer, generating an overall stronger effect compared 
to the same force evenly distributed on the monomers. 
In fact, a force generates a stronger effect if it is placed 
in the region near the free polymer end (its effect is re- 
distributed also to all preceding bonds). The curves in 
Fig|6]and Figl7]have been obtained with the theoretical 
formulations presented in this Section and confirmed by 
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TABLE I. Asymptotic forms of the force-extension curves for 
all cases described in the paper: FJC and WLC models in 2D 
and 3D geometry with force applied / or field applied g. 



Asymptotic form 
Polymer chain of for /, g — > 



Equation 



Jd- or ^ 

knT k„T 



Asymptotic form 
of 7^ for f,g ^ CO 



T — or 

kpiT kp,T 



regime) and 1 < log(7V + 1) < N (high force regime) for 
the FJC model and 1 < 1 + N/2 < N (low force regime) 
and 



/iV < 2{y/N + 1 -1) < N (high force regime) for 
the WLC model (always for > 2). 



V. ACTION OF A PULLING FORCE NOT 
ALIGNED WITH THE EXTERNAL FIELD 



FJC (2D) / 

Eg. QUI 

FJC (3D) / 
Eo.i3T\ 




FJC (3D) g 
Bo.l33l 

WLC (2D) / 

^ J -' 

WLC (3D) / 
WLC (2D) g 
WLC (3D) g 



1 



1 

— X 

3 



X 



l + -^x 



1 + 



log(jV + 1) 1 
27V X 



log(Af + 1) 1 
N X 



1 



1 1 



4 /I 



2 L/p 
3T" 



1- 



1 1 



2 /I 



1 + 



N 



2 Z/p 
3T 



1 + 



iV 



X 1 



X 1 



I 



1 /?vTT- 1 



27V 



1 yWTT-l 



N 



a series of MC simulations. In all case we obtained a 
quite perfect agreement between the two formulations. 
The knowledge of the closed-form expressions allowed us 
to analytically analyze the behavior of the chains for very 
low and very high applied forces (or fields). The results 
are shown in Table |TJ interestingly, we note that the ex- 
tension is always a linear function of the small applied 
perturbation. Nevertheless, the corresponding constant 
of proportionality depends on TV only when a field is ap- 
plied to the chain; conversely, it is independent of TV with 
a single force applied at one end. On the other hand, with 
a large perturbation applied to the molecule, we observe 
a 1/a; behavior for the FJC models and a l/i/a: behavior 
for the WLC models. To conclude we also remark that 
the order of the curves observed in Fig Eland Fig d is con- 
firmed also in the low and high force (or field) regime by 
the following inequalities: 1 < 1 -I- N/2 < N (low force 



In previous Sections we considered the polymer chain 
immersed in an external field with an external force equal 
to zero at its end. However, since we developed a form of 
the partition function also taking into account an exter- 
nal force applied at the end of the chain (at least for the 
FJC model), we can directly study the important case 
with a non zero force superimposed to an external field, 
in general having different orientation. To do this, we 
keep fixed the origin of the chain and apply a constant 
force at the end of the polymer with different angles with 
respect to the direction of the applied field. We will anal- 
yse such a problem for both the FJC and WLC cases. 

To begin, we consider a pulling force perpendicular to 
the direction of the applied field, respectively the y and z 
axis of our reference frame. For increasing values of the 
bending spring constant k going from nearly zero (FJC 
model) to 8 • 10~^^Nm (WLC model, including the bend- 
ing constant of the DNA given by k = 0.6 • 10~^^ Nm 
~ ISfcfiT). In FiglHwe reported the results for the av- 
erage monomers positions and their variances. The red 
solid lines correspond to the analytical results for the 
FJC case, while the black symbols correspond to the MC 
simulations. It is interesting to observe the effect of the 
persistence length (or, equivalently of the bending stiff- 
ness): in fact, in the top panel of FigjS] we note that 
the chains with an higher bending spring constant tend 
to remain more straight under the same applied load. 
At the same time, in the fourth panel of FigjS] we ob- 
serve a decreasing variance along the z-axis (direction of 
the applied field) with an increasing bending spring con- 
stant; this fact can be easily interpreted observing that 
an higher rigidity of the chain reduces the statistical fluc- 
tuations in the direction of the applied field. The situ- 
ation is more complicated for the variances along the x 
and y directions: in fact, along the chain, there are some 
monomers with variances larger than the corresponding 
FJC case and others with smaller values. 

In Fig 12] the average positions of the monomers for dif- 
ferent directions of the external force are reported. The 
figure shows how the average monomer positions depend 
on the bending rigidity k and on the external force angle 
6. As before we can observe that the persistence length of 
the chain tends to maintain a low curvature in the shape 
of the chain. This phenomenon is more evident with an 
increasing angle between the force and the field. In fact, 
in FigEl the deviation between the FJC results and the 
WLC ones is higher for the angles approaching tt, where 
the force and the field are applied in opposite directions. 

In FigsdU] and [11] the three components of the vari- 
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FIG. 8. (color online) Action of a pulling force / (along 
the y-axis) perpendicular to the applied field g (along the 
z-axis). We adopted different values of the bending spring 
constant: k = 0.08,0.6,2,8 ■ 10"^^ Nm. The chain length is 
fixed [N = 20), the external field amplitude is g = 4 pN and 
the force applied to the last monomer of the chain corresponds 
to / = 8 pN. The red solid lines correspond to the analytical 
results for the FJC case (see Eas. ((29)) and (|42)) '). Black circles 
correspond to the MC simulations with the different bending 
spring constants. In the top panel we reported the average 
positions, while in the others the three variances of the x, y 
and z components. 



FIG. 9. (color online) Average positions of the chain for dif- 
ferent angles between the external traction force / and the 
direction of the applied field g. We adopted A'^ = 20, g = 4 pN 
and / = 60 pN. The red solid lines correspond to the FJC ana- 
lytical result, Ea. (|29p . The symbols represent the MC results 
for the WLC model with k = 0.08, 0.6, 2 ■ 10"^^ Nm (circles, 
triangles and squares, respectively) . For both FJC and WLC 
models we used different values of the angle between the ap- 
plied field and the traction force 9 = 7r/2, 37r/4, 57r/6, 157r/16 
from the right left. 



ance are reported versus the position of the monomer 
along the chain and the angle between the field and the 
force directions, for the FJC and WLC case, respectively. 
We can extract some general rules about this very com- 
plex scenario: as for the variance along the x direction 
we observe it to be an increasing function both of the po- 
sition i along the chain an of the angle between / and 
g. Both behaviors can be interpreted with the concept 
of persistence length, as discussed above. Conversely, 
the description of the variance along the y direction is 
more complicated. In fact, while the increasing trend of 
the variance with the position i along the chain is main- 
tained, we observe a non monotonic behavior in terms of 
the angle 0, with a minimum of the variance at about 
9 = 2tt/3. Finally, the variance along the z direction is 
always increasing along the chain, but it shows a maxi- 
mum near 9 = t: (at least in the first part of the polymer 



chain) . 



VI. CONCLUSIONS 

In this work we investigated mechanical and confor- 
mational properties of fiexible and semi-fiexible polymer 
chains in external fields. As for the FJC model we de- 
veloped a statistical theory, based on the exact analyt- 
ical determination of the partition function, which gen- 
eralizes previous results to the case where an external 
field is applied to the system. In particular we obtained 
closed form expression for both the average conformation 
of the chain and its covariance distribution. For sake of 
completeness, all calculations have been performed both 
in two-dimensional and three-dimensional geometry. On 
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i/N 00 9 



i/N 9 



i/N 9 



FIG. 10. (color online) Monomer variances versus the position along the chain (i) and the angle between force and field 
(0 < 6* < tt) for the FJC model. As before we used iV = 20, p = 4 pN and / = 60 pN. 




i/N 9 i/N 9 i/N o o 9 

FIG. 11. (color online) Monomer variances versus the position along the chain (i) and the angle between force and field 
(0 < ^ < tt) for the WLC model. As before we used A'' = 20, g ~ 4 pN and / = 60 pN. We also adopted a bending stiffness 
K = 0.6 • 10"^*' Nm. 



the other hand, as for the WLC model we derived new 
approximate expressions describing the force-extension 
curve under the effect of an external field. They can be 
considered as the extensions of the classical Marko-Siggia 
relationships describing the polymer pulled by a single 
external force applied at the free end of the chain. All 
our analytical results, for both FJC and WLC models, 
have been confirmed by a series of Monte Carlo simu- 
lations, always found in very good agreement with the 
theory. 

The overall effects generated on the tethered polymer 
by the application of an external field can be summarized 
as follows. As for the average configuration of a chain, 
it is well known that a single pulling force generates a 
uniform deformation along the chain (for a homogeneous 
polymer with all monomers described by the same effec- 
tive elastic stiffness). On the contrary, the application 
of an external field produces a non uniform deformation 
along the chain, showing a larger deformation in the por- 
tion of the chain closest to the fixed end. Moreover, the 
variances of the positions increase linearly along the chain 
with a single force applied to the polymer. Conversely, 
the polymer subjected to an external field exhibits a non- 



linearly increasing behavior of the variances along the 
chain. More specifically the variances assume the largest 
values nearby the last free monomers, where we can mea- 
sure the highest fluctuations. 

To conclude, we underline that the use of the MC 
method, once validated against known analytical solu- 
tions, is crucial for analysing models conditions which 
are beyond reach of a full analytical calculation. We take 
full profit of this approach for analysing the effects of the 
combination of an applied force at the free end together 
with an external field, especially when the two are not 
aligned. We have analysed the average configurational 
properties of the polymer, observing a very complex sce- 
nario concerning the behavior of the variances. 
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